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The species doubling problem of the lattice fermion is resolved by introducing hopping interactions that mix 
left- and right-handed fermions around the momentum boundary. Approximate chiral symmetry is realized on 
the lattice. The deviation of the fermion propagator from the continuum one is small. 



1. INTRODUCTION 

In contrast with the great success of lattice 
gauge theory, lattice fermions remain a long- 
standing problem. Naive discretization causes the 
species doubling problem pQ. The situation does 
not change regardless of how the lattice spacing 
is reduced as long as the space-time derivative is 
modeled as a naive difference. Many attempts 
have been made to fix the doubling problem 
1 1121,'il ll:")l()l7j . Wilson removed doublers at low 
energy by introducing an interaction that mixes 
left- and right-handed fermions 1 . However, un- 
wanted degeneracy persist at high energy and chi- 
ral symmetry is explicitly broken. To fix these 
problems, Kaplan modified Wilson's fermion with 
an extra dimension and succeeded in realiz- 
ing approximate chiral symmetry. However, the 
cost of numerical calculations based on it is not 
cheap. If we find a method to perform such calcu- 
lations without the extra dimension, calculation 
time decreases largely and a deeper understand- 
ing of quantum field theory becomes possible. 

The lattice fermion has another serious prob- 
lem. The fermion propagator defined on a lat- 
tice deviates from the continuum one even if the 
doublers are removed with the existing techniques 
such as Kaplan's fermion [5]. We need to modify 
the discretized propagator somehow so that it is 
close to the continuum one as far as possible. 

In general, the extra dimension can be ex- 
pressed as hopping interactions in a lower- 
dimensional system. Kaplan's fermion is a formu- 
lation with an extra dimension, so there must be a 
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corresponding Hamiltonian with no extra dimen- 
sion. Also, the shape of the fermion propagator 
can be improved with hopping interactions. 

In this talk, based on a Hamiltonian formal- 
ism, we introduce ultralocal hopping interactions 
to remove doublers and improve momentum de- 
pendence of fermion energy. (The word "ultralo- 
cal" means that fermion hopping is restricted to 
a finite range on a real-space lattice From 
knowledge of the continuum theory, we know the 
correct momentum dependence of the energy. We 
start from momentum space and go back to real 
space by way of discrete Fourier transform. See 
Ref. P| for the details of this work. 



2. FORMULATION 

Let us consider a free Dirac fermion on a (1+1)- 
dimensional Hamiltonian lattice. 

N/2 

H=- J2 wCitVi, (1) 

l=-N/2+l 

where I is an index for momentum and N is the 
number of sites. Q an d are discrete Fourier 
transform of real-space two-component fermion 
operators ifj n and ip n = ip^j for n = 1, 2, . . . , N, 
respectively. The system is quantized with anti- 
commutators in the usual way. Periodic bound- 
ary conditions are assumed in real space, pi = 
2ttI/N. We try to create a real-space Hamiltonian 
with no doubler that reproduces Eq. . To find 
real-space representation of pi , let us consider the 
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following function 



a=l 



a 



sinap. 
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In the limit M — ► oo, the function @ goes to 
p. The function necessarily has a node at the 
boundary p = ±7r because it has a periodicity of 
2ir. The node is the cause of the doubling prob- 
lem. The doubler remains as a singularity at the 
boundary even if the limit M — > oo is taken (by 
"doubler" we mean unwanted energy degeneracy 
that is not contained in the continuum theory). 
Anyway, the parameter M needs to be small for 
practical formulation because M corresponds to 
the maximum distance of fermion hopping. 
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Figure 1. The solid line plots the correct energy 
±p from the continuum theory. The dashed line 
plots the function ±s(p) for M — 5 with the Lanc- 
zos factor. The dot-dashed line plots the function 
±k(p) for M = 5, u = 130, and v = 8.4 with the 
Lanczos factor. The function ±k(p) almost agrees 
with ±p in the fundamental region \p\ < tt except 
for a small deviation around momentum \p\ ~ 2.3. 



In addition to the doubler around the momen- 
tum boundary, the function (J2J) has another de- 
generacy. It oscillates around p and has local min- 
ima if the summation is truncated with a small 
M. In Fourier analysis, it is called the Gibbs 
phenomenon, which occurs if a function to be ex- 
panded has a singularity ^Oj ■ The oscillation can 
be removed by replacing Eq. (J2J) with 
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s (p) = E S <* sin ap, 

a=l 

where 
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F a is called the Lanczos factor |10|. The factor al- 
most removes the oscillation of Eq. (J2J. However, 
the doubler modes around the boundary still re- 
main. We are going to remove them by a trick 
with hopping interactions. Let us consider the 
following momentum-space Hamiltonian: 
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l=-N/2+l 
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where s/ = s(pi)/a and m is fermion mass. We 
introduce interactions cj that mix left- and right- 
handed fermions 
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where q are assumed to be nonzero only for 
N/2. The Hamiltonian © can be diagonalized. 
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where ki = y/sf + (m + q) 2 are energies of one 
particle states and £/ are transformed variables. 
The coefficients q are determined so that energy 
of one-particle states becomes close to the contin- 
uum one as far as possible. 

In Fig. ^ the functions ±p, ±s(p), and ±k(p) 
are compared, where k{p) = y/ s 2 (p) + c 2 (p). The 
function ik(p) corresponds to energy of one- 
particle states given by the Hamiltonian JSJ) with 
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to = 0. The function ±k(p) agrees well with the 
correct energy ±p from the continuum theory in 
the fundamental region |p| < w. The Hamiltonian 
|5| approximately reproduces the continuum the- 
ory without doubler if we identify c; = c(pi)/a 
and ki = k(pi)/a. 
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Figure 2. The matrix elements of 75 are plot- 
ted as functions of the momentum pi for a = 1, 
M = 5, u = 130, and v = 8.4 with the Lanczos 
factor. For a lattice size N = 40, the functions 
are plotted in the right-half plane of momentum 
space pi > 0. The circles are plots of the diagonal 
(1,1) element of 75. The crosses are plots of minus 
of the off-diagonal (1,2) element of 75. The solid 
and dashed lines are plotted to guide the eyes. 



Figure shows the diagonal and off-diagonal 
matrix elements of the transformed 75 in the new 
basis that diagonalizes the Hamiltonian (JSJ) with 
m = 0. The diagonal (1,1) element of 75 is almost 
unity at low and intermediate energy pi < 2.3 
and deviates from unity at p\ > 2.3. The off- 
diagonal (1,2) element of 75 oscillates around zero 
at p\ < 2.3 and becomes unity at pi = tt. At low 
energy, the deviation of the off-diagonal elements 
from zero is not large and becomes smaller as the 



parameter M increases. At p\ < 2.3, the trans- 
formed left- and right-handed fermions have ap- 
proximately the correct chiral charges 1 and — 1, 
respectively. The low-energy Hamiltonian (J5J has 
approximate chiral symmetry because the com- 
mutation relation between the Hamiltonian JSJ 
and chiral charge defined with 75 is almost zero 
for small The errors associated with chiral sym- 
metry can be improved in a systematic way by 
increasing M. The value used here for the pa- 
rameter M is sufficiently small and does not deny 
application of the model to actual numerical anal- 
ysis with a computer. 

The real-space Hamiltonian and the Euclidean 
action for the Dirac spinor is obtained by substi- 
tuting the discrete Fourier transform into Eq. J5J • 
As usual, gauge symmetry can be implemented by 
inserting exponentiated gauge fields between the 
hopping fermions. The continuum limit a — > is 
taken with M fixed. The method can be extended 
to higher-dimensional cases by taking care of the 
doubler of each direction in the same way. 

In this talk, I have shown how to realize ap- 
proximate chiral symmetry at low energy. Ex- 
plicit breaking of chiral symmetry has been com- 
pressed to high energy. In future works, it should 
be precisely checked if insertion of gauge interac- 
tions affects chiral properties. 
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